On the appearance of hyperons in neutron stars 



O 

oo 



> 
cn 

On 



00 

o 

> 

X 
J3 



H. Dapo,i'B B.-J. Schaefer,2.[t] and J. Wambachi'^ 

^Institut fiir Kernphysik, TU Darmstadt, Schlofigartenstr. 9, D-64289 Darmstadt, Germany 
^ Institut fur Physik, Karl-Franzens-Universitdt Graz, Universitdtsplatz 5, A-8010 Graz, Austria 
^ Gesellschaft fiir Schwerionenforschung mbH, Planckstr. 1, D-64291 Darmstadt, Germany 

(Dated: February 8, 2010) 

By employing a recently constructed hyperon-nucleon potential the equation of state of /3- 
equilibrated and charge neutral nucleonic matter is calculated. The hyperon-nucleon potential is 
a low-momentum potential which is obtained within a renormalization group framework. Based 
on the Hartree-Fock approximation at zero temperature the densities at which hyperons appear in 
neutron stars are estimated. For several different bare hyperon-nucleon potentials and a wide range 
of nuclear matter parameters it is found that hyperons in neutron stars are always present. These 
findings have profound consequences for the mass and radius of neutron stars. 



PACS numbers: 13.75.Ev, 21.65. Mn, 26.60.-c 

I. INTRODUCTION 

Neutron stars (NS) are compact objects with interior 
densities of several times normal nuclear density. The 
precise and detailed structure and composition of the in- 
ner core of a NS is not known at present. Several possibil- 
ities such as mixed phases of quark and nuclear matter, 
kaon or pion condensates or color superconducting quark 
phases are under debate. 

In the present work the influence of hyperons with 
strangeness S = —1 (A, E~. S° and E"*") on the composi- 
tion and structure of a NS is investigated. In this context 
the central and essential quantity to be analyzed is the 
equation of state (EoS). The EoS determines various NS 
observables such as the mass range or the mass-radius 
relation of the star. The composition and structure of 
neutron and also of protoneutron stars have been investi- 
gated in detail with a wide range of EoS for dense nuclear 
matter [l|, 01 . The emergence of hyperons for increasing 
nucleon densities has been suggested in the pioneering 
work of Q . Since then the impact of hyperons on dense 
matter has been studied extensively with different ap- 
proaches, see e.g. P, 0-11] Unfortunately, the details of 
the hyperon-nucleon (YN) interaction and even more of 
the hyperon-hyperon (YY) interaction are known only 
poorly. The limited amount of available experimental 
data enables the construction of many different poten- 
tials. For example, the Nijmegen group has proposed 
six different potentials which all describe the low-energy 
data such as phase shifts equally well, see e.g. (To| . 

We wish to explore the differences between the avail- 
able YN interactions and their influences on the appear- 
ance of hyperons in neutron stars. For this purpose the 
hyperon YN threshold densities are calculated. Since the 
5]° and E+ hyperons are heavier than the A and E~ hy- 
perons they typically appear later. Hence, we will focus 
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on the A and E~ hyperons in the following. 

The paper is organized as follows: In Sec. lIIII the EoS of 
dense matter including hyperons are calculated. For the 
pure nucleonic part of the EoS a parameterization is used. 
A central quantity which enters the EoS is the single- 
particle potential. The derivation of the single-particle 
potential for hyperons is given in Sec. ]llj The require- 
ments which are necessary for equilibrium are discussed 
in Sec. IIV Al and the results are collected in Sec. IIVBI 
The next section Sec. ]V] shows the calculations of neu- 
tron stars with and without hyperons. Finally, the work 
is summarized and conclusions are drawn in Sec. IVII 



II. SINGLE-PARTICLE POTENTIALS 

Single-particle potentials are a useful and important 
tool to determine the density at which hyperons begin to 
appear in baryonic matter. They can be obtained from 
an effective low-momentum YN potential T^ow k [HI, 
in the Hartree-Fock approximation. Details concerning 
the derivation of the YN single-particle potentials in this 
approximation can be found in [l^. With these single- 
particle potentials the chemical potentials and particle 
energies in /3-equilibrated matter can be calculated. This 
allows to compute the threshold densities for the appear- 
ance of a given hyperon species and to establish the con- 
centrations of all particles in dense matter at a given den- 
sity. In the process we will also determine the EoS for 
the mixture of leptons and baryons (electrostatic inter- 
actions are neglected because their energies are orders of 
magnitude smaller than the other interaction energies). 

Our many-body scheme employs several "bare" YN in- 
teractions as input for the Viow k calculation: the origi- 
nal Nijmegen soft core model NSC89 1141 . a series of new 
soft core Nijmegen models NSC97a-f [lOll . a recent model 
proposed by the Jiilich group J04 [T^ and chiral effec- 
tive fleld theory (xEFT) [l^. The first three models are 
formulated in the conventional meson-exchange (OBE) 
framework, while the xEFT is based on chiral pertur- 
bation theory (for recent reviews see e.g. |17Hl9f ). In 
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the xEFT approach a cutoff A^eft enters which we fix 
to A^EFT = 600 MeV and label the results obtained by 
XEFT600. 

In general, the effective low-momentum Y N interac- 
tion, the V\o^N k for hyperons, is obtained by solving 
a renormalization group equation. The starting point 
for the construction of the Viow k is the half-on-shell 
T-matrix. An effective low-momentum Tjow k-nratrix 
is then obtained from a non-relativistic Lippmann- 
Schwinger equation in momentum space by introducing 
a momentum cutoff A in the kernel. Simultaneously, the 
bare potential is replaced with the corresponding low- 
momentum potential Viow k, 
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The on-shell energy is denoted by while q', q are rela- 
tive momenta between a hyperon and nucleon. The labels 
y, J/', z indicate the particle channels, and a, a', /3 denote 
the partial waves, e.g. a ~ LSJ where L is the angular 
momentum, J the total momentum and S the spin. In 
Eq. ([1]) the energies in the denominator are given by 
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MyMN/My, where My is 
the total mass of the hyperon-nucleon system. My = 
My + Mn. Finally, the effective low- momentum Viow k 
is defined by the requirement that the T-matrices are 
equivalent for all momenta below the cutoff A. Details 
for several bare YN potentials within the RG framework 
can be found in Refs. [Til [T2|. 

From the different effective V|ow k interactions we cal- 
culate the single-particle potential Ub{p) of a baryon 
b G {p, n. A, S", S", E!~, 5°} with three-momentum 
p ^ \p\. In general, it is defined as the diagonal part 
in spin and flavor space of the proper self-energy for 
the corresponding single-particle Green's function. In 
the Hartree-Fock (HF) approximation for a uniform sys- 
tem it represents the first-order interaction energy of the 
baryon with the filled Fermi sea. It is evaluated as the di- 
agonal elements of the low-momentum potential matrix, 
V^{q)^ where an evident short- hand labeling for the di- 
agonal elements has been introduced (cf. Eq. ([1])). Note, 
that the relative momentum is given hy q ~ \p — p'\ and 
the particle channel index is given by y = bb' . In the 
HF approximation the single-particle potential has two 
contributions: the (direct) Hartree- and the (exchange) 
Fock-term f20| 
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In this expression the diagonal elements of the nucleon- 
nucleon (NN) interaction is also included. In principle. 
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FIG. 1: Density dependence of A single-particle potentials 
for various hyperon-nucleon interactions in symmetric nuclear 
matter. 



an effective low-momentum potential for the NN inter- 
action is also known but we will use for the A^A^-sector a 
parametric Ansatz to be discussed later (see Eq. ©). 
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FIG. 2: Same as in Fig. [1] but for the E~ hyperon. 



The density dependence for several A potentials at rest 
in symmetric nuclear matter (no hyperons present) is 
shown in Fig. [TJ The square represents the generally ex- 
cepted empirical depth of Ua{p = 0) « —30 MeV. While 
most of the potentials used can reproduce this value, 
the Jiilich potential (J04) yields a stronger binding while 
the old Nijmegen potential (NSC89) underestimates the 
binding. All other potentials agree up to the saturation 
density. However, with increasing density, the differences 
grow, leading to different bindings at rest. This will have 
consequences for the predictions of the A hyperon con- 
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centration in dense nuclear matter. A comparison with 
other works, [2l|, and [T3| shows some differences 
but mostly yields similar results. 

Fig. [5] shows the density dependence of several po- 
tentials at rest in symmetric nuclear matter similar to 
Fig. [H No agreement of the various U-^- potentials over 
the density range considered is seen. Compared with a G- 
matrix calculation a stronger binding for the S" single- 
particle potential is obtained. In Ref. [l^l further details 
concerning the hyperon single-particle potentials can be 
found. 

For constructing the total energy /particle one has to 
distinguish which part of the total single-particle poten- 
tial comes from the in-medium nucleon interaction and 
which from the hyperons. For this purpose, the single- 
particle potential can be split into a nucleonic and hy- 
peronic contribution 



U,{p)^U,^{p) + U,^ip) 



(4) 



where the first one, U^(ji), denotes the baryon interact- 
ing with a nucleon and the second one, {p), means the 
baryon interaction with a hyperon. Note, that the baryon 
can either be a nucleon or a hyperon. Accordingly, the 
nucleonic contribution is calculated from Eq. ([3]) via 
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and analogously for the hyperonic contribution. 



III. EQUATION OF STATE 

By means of the single-particle potential the total en- 
ergy per particle, E/A, can be easily calculated. It is 
given by the total energy of a baryon of mass Mi, and its 
kinetic and potential energy Uh(j>) divided by the total 
baryon number density pB'- 
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The potential Ub{p) describes the average field which 
acts on these baryons due to their interaction with the 
medium. The baryon Fermi momentum kp^ is given by 



(7) 



with the baryon fraction ratio Xb = Pb/PB for baryon b. 

It is well-known that non-relativistic many-body cal- 
culations, based on purely two-body forces, fail to re- 
produce the properties of nuclear matter at saturation 
density. This is also the case in the present work. In 
order to proceed we replace the purely nucleonic contri- 
butions (without the influence of the hyperons) by an 
analytic parameterization developed by Heiselberg and 



Hjort- Jensen \2^. This replacement makes the study of 
hyperons more robust since the NN sector can be con- 
trolled more easily. 

From Eq. ^ the purely nucleonic contribution to the 
energy per particle reads 
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where the single-particle potential U^{p) only contains 
the nucleonic contribution, cf. Eq. ([S]). However, instead 
of using the Hartree-Fock expression we employ the fol- 
lowing parameterization 



Enn/An = Mn - Equ- 
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where u ~ Pn/po denotes the ratio of the total nucle- 
onic density pn = {xp + Xn)pB to the nuclear satura- 
tion density po = 0.16 fm"'^. The corresponding proton 
and neutron fraction are denoted by Xp and a;„. The 
parameters Eq , S, Sq are related to properties of nuclear 
matter at saturation density, i.e. Eq is the binding en- 
ergy per nucleon at saturation density while Sq and 6 are 
connected to the symmetry energy and incompressibility, 
respectively. 

As mentioned in the previous section, cf. Eq. ([3]), we 
can separated the potential contribution of nucleons into 
one coming from the interaction with other nucleons 
U^{p), and one coming from the interaction with hyper- 
ons Ujj{p). The latter does not contribute to the purely 
nucleonic EoS when we consider symmetric matter only. 
However, in the next section, when we investigate hyper- 
ons such terms are considered. 

The parameterization ^ is fitted to the energy per 
particle in symmetric matter obtained from variational 
calculations with the Argonne Vig, nucleon-nucleon inter- 
action including three-body forces and relativistic boost 
corrections [2j|. The best fit parameters are Eq = 
-15.8 MeV, 50 = 32 MeV, 7 = 0.6 and 5 = 0.2. The 
EoS from Ref. [IJl is considered as one of the most reli- 
able ones. In this way possible uncertainties coming from 
the nucleonic EoS are minimized. 

The symmetry energy in dense matter is defined as 



1 d'^E/A 
8 dxl 



(10) 
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Since there are no hyperons at saturation density we can 
use Eq. (jH]) directly to obtain at — Sq- The incompress- 
ibility is given by 



9p- 



d^E/A 



(11) 



from which we obtain the relation Kq = — 18i?o/(l + ^). 

To study the impact of the softness of the EoS and 
the effects of the symmetry energy we vary Kq and at 
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in a broader interval. From experimental constraints the 
values for Kq range between 200 MeV and 300 MeV and 
those for at between 28 MeV and 36 MeV, see Ref. ^ 
and references therein. We vary Kq and at within these 
limits to study the effects of both parameters on the ap- 
pearance and concentrations of hyperons in dense matter. 
While they do not influence the particle concentrations 
directly they modify the composition of the matter indi- 
rectly by changing the available energy. In this way the 
point at which hyperons will appear is affected. Results 
for various values of Kq are shown in Fig. [3] (in symmet- 
ric matter the energy per particle is only sensitive to the 
incompressibility) . In all cases the saturation point is at 
E /A ^ —16 MeV. The parameter range allows us to clas- 
sify the nucleonic EoS as a stifFer {Kq = 300 MeV) or a 
softer (A'o = 200 MeV) one. 



with the rest 




In symmetric matter, which is composed only of nu- 
cleons, E' /A vanishes, but with this separation we can 
calculate the total E/A for arbitrary liyperon concen- 
trations. In the following we will calculate the EoS in- 
cluding hyperons by determining their concentrations in 
/3-equilibrium. 
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FIG. 3: Parametric EoS in symmetric nuclear matter as ob- 
tained from the parameterization of Ref. [2^. 



IV. /J-EQUILIBRIUM 

A. Composition of matter 

In order to determine the threshold densities for hy- 
perons their concentrations are needed. These are fixed 
by charge neutrality and /3-equilibrium. The latter refers 
to the equilibrium under the weak interaction decays 



Bi^ B2 + I + D1 



(14) 



where Bi and B2 denote the baryons, I e {e^,/i~,r~} 
the negatively charged leptons and vi the corresponding 
neutrinos. In the case when the neutrinos are not trapped 
in the star (i.e. /i^ = 0) these requirements amount to 
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In addition, Kq directly influences the maximum al- 
lowed mass of a neutron star. By increasing Kq the en- 
ergy of the system is increased and as a consequence, 
more and more hyperons can be produced. This in turn 
will decrease the allowed maximum mass of a neutron 
star. Such a nontrivial connection creates a conundrum: 
if we use a stifFer nucleonic EoS by increasing Kq we then 
allow for higher hyperon concentrations which softens the 
total EoS. 

By means of Eq. ([8|) we can split the total energy per 
particle, Eq. ([S]), into a purely nucleonic part and a re- 
mainder E' /A via 
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for the charge neutrality and 



/iE- = = ^^n+ fJ-e, 

Mho = = = Pn, 

/^S+ = Pp = Pn — Me, 



(16) 
(17) 
(18) 



for the chemical potentials. The densities of positively 
and negatively charged baryons and leptons are denoted 
by and p[^\ respectively. The chemical potentials fi 
are labeled by the corresponding particles. In the absence 
of neutrinos all lepton and antilepton chemical potentials 
are equal. In addition to the electrons, muons are also 
present. The r-lepton does not appear since it is too 
heavy. 

At zero temperature, the chemical potential of a 
fermion system is equal to its Fermi energy. For rela- 
tivistic noninteracting leptons it is given by 



Pi = 



(19) 
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FIG. 4; Density ratios for different particles for a "soft" nu- 
cleonic EoS as a function of the baryon density using the 
XEFT600 model. 



FIG. 5: Same as Fig.|4]but for "stiff" nucleonic EoS using the 
NSC97f model. 



with the corresponding lepton density ratio xi = pi/pL- 
The total lepton density is the sum over all three lep- 
tons. For nonrelativistic interacting baryons, the chemi- 
cal potential for species b reads 



fJ-b 
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For a given total baryon density ps the equations (|15p- 
()18p govern the composition of matter, i.e. the baryonic 
and leptonic concentrations. The corresponding solution 
is referred to as /3-stable matter. 

For the sake of consistency we now have to treat the 
nucleonic part of the chemical potential /ijy in the same 
way as the corresponding energy per particle. Since the 
chemical potential can be obtained as a derivative of the 
energy density e and is related to the energy per particle 
via e ~ pbE/A, we use the definition 



Pb = 



dpb 
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to have the appropriate replacement in the nucleonic 
chemical potential. Finally, we arrive at the expression 
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where we have effectively replaced Mjv + jm]^ + ^Ni^F^ ) 
of Eq. with the derivative deNN/dpN- In this way 
the parameterization Eq. ([9]) enters in the nucleonic part 
of the chemical potential. 

Since we are only parameterizing the nucleonic sec- 
tor, no such replacement is necessary for the hyperons. 
However, since we have neglected the YY interaction 
Uy (kFY) is zero and Eq. ([20| reduces to 
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For the determination of the particle concentration the 
single-particle potential in equilibrium is used. For hy- 
perons below the threshold density it is given by Uyip = 
0), similar to the symmetric matter case. Above the 
threshold density it depends on composition and den- 
sity. In Fig. [6] the density dependence of U-^- (kp^ ) in 
/3-equilibrium for two different incompressibilities Ka is 
shown. In the figure a kink in the curves appears at the 
point where the hyperons appear. 

Another observation is the relative ordering and the 
magnitudes which resemble those of the single-particle 
potentials at zero momentum in symmetric matter as 
shown in Fig. [H Essentially, the NSC97a, NCS97c, 
NSC97f and J04 interactions are still slightly attractive 
while the NSC89 and xEFTGOO remain repulsive. Simi- 
lar observations hold for the A system. A new structure 
in form of a second inflection point emerges due to the 
appearance of the hyperon. 

A better indicator at which densities hyperons start to 
appear is given by the concentrations of all particles and 
is displayed in Figs. 2] and [5] In Fig. U] a "soft" nucle- 
onic EoS is used in combination with an attractive AA^ 
and a very repulsive T,N interaction implemented by the 
XEFT600 model. In contrast in Fig. [5] a "stiff" EoS is 
used represented by the NSC97f model which has a simi- 
lar AA^ interaction compared to the xEFT600 model but 
also an attractive J^N interaction. This difference already 
leads to very different density profiles. While in Fig. |4] 
the A hyperon is the first one which appears and no 
hyperons are present, the hyperon appears first in 
Fig. [51 

One should note that with the appearance of the 
hyperon the density of the negatively charged leptons 
starts to drop immediately. This is because their role in 
the charge neutrality condition, Eq. (jlSp . is now being 
taken over by the E~. Similarly, the appearance of the 
A hyperon will accelerate the disappearance of neutrons 
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FIG. 6; Density dependence of U^- {kF^_ ) for /3-eqmlibrated 
matter. Upper panel: Kq = 200 MeV, lower panel: Kq = 300 
MeV. 




FIG. 7: Various EoS for /3-equilibrated matter for two differ- 
ent Ko as a function of density. Upper panel: soft EoS, lower 
panel: stiff EoS. 



since both are neutral particles. 

Once the composition of the matter has been deter- 
mined by demanding /3-equilibrium we can calculate the 
energy per particle. For this purpose, we cannot use 
Eq. dl]), but have to use Eqs. and The re- 

sult is presented in Fig. [7] where the energy per particle 
in /^-stable matter is shown as a function of the den- 
sity for different YN models. The symmetry energy is 
fixed to at = 32 MeV while the incompressibility is set 
to Kq ~ 200 MeV (upper panel in the figure) and to 
Kq = 300 MeV (lower panel). In addition, the EoS with 
hyperons is compared with the purely nucleonic one. 

One easily observes the onset of the hyperon appear- 
ance at the point at which the curves start to deviate. As 
expected the differences between the various YN interac- 
tions do not modify the EoS for very small densities. In 
the range between (2 — 3)po, all EoS's are similar to each 
other. However, for increasing densities the influence of 
hyperons becomes more significant resulting in rather dif- 
ferent EoS's. This concerns not only the magnitudes of 
the different energies per particle but also their slopes at 
higher densities. These variations will lead to differences 
in the pressure and finally to significant changes in the 



possible maximum mass of a neutron star. 



B. Threshold densities 

The appearance of a given hyperon species is deter- 
mined by increasing the density for fixed Kq and at . The 
resulting threshold densities for the S" hyperon for cer- 
tain Kq and at are collected in Fig.[5]for six different YN 
interactions. Similarly, the threshold densities for the A 
hyperon are shown in Fig. [5] From these figures one sees 
how the single-particle potentials for various YN inter- 
actions modify the threshold densities. In this way, the 
properties of the YN interaction in Fig. [T] and Fig. [2] can 
be attributed to the hyperon appearances. 

From Fig. [8] one concludes that the E" hyperon ap- 
pears between 1.4po and 2Apo with the exception of the 
xEFTGOO model. For almost all YN interactions used in 
the present study the E~ is the first hyperon which will 
appear even though the A hyperon is the lighter one. The 
reason is that the heavier mass of the E" is offset due to 
the presence of the e~ chemical potential, cf. Eq. ([T5|) . 
In general, heavier and more positively charged particles 
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FIG. 8: Threshold densities for the E~ hyperon with various incompressibilities Kq and symmetry energy at for six different 
YN interactions. 



appear later. In the case of the S^, compared to the 
A, the effect caused by the electric charge dominates the 
one coming from the mass in almost all cases. For the 
hyperon a further modification caused by the electric 
charge, is the influence of on the threshold density be- 
cause the electron chemical potential is modified by the 
symmetry energy. Thus, the decrease of the threshold 
densities due to the increase of Kq is analogous to the 
increase due to at- 

For the A hyperon the range of threshold densities is 
between densities from 1.7po to 4.5po depending on the 
choice of Kq, at and the YN interaction used, cf. Fig. [HI 
The influence of Kq on the threshold density for this hy- 
peron is larger than the one from at- This is reasonable 
since Kq controls the rate of the energy increase with the 
density more directly, while at affects only the details of 
the /3-equilibrium. One clearly recognizes in Fig. [9] that 
the A appears earlier for larger incompressibilities. Thus, 
in general we see that for increasing A'o the threshold 
densities decrease for both hyperons. 

In contrast to the previous and at discussion, the 



influence of the single-particle potentials on the threshold 
densities is harder to analyze. 

The threshold densities for the E~ are largest for the 
XEFT600 interaction which yields the most repulsive 
single-particle potential. In general, hyperons will ap- 
pear earlier for a more attractive single-particle poten- 
tial. This becomes obvious from Eq. (|20p : the chemical 
potential decreases for a more negative Ub{kp^) and, con- 
sequently, the threshold density will also decrease. Thus, 
the most repulsive single-particle potential like the one 
for the xEFT600 leads to the largest threshold density. 
For the A hyperon the threshold densities are smallest 
for the most attractive single-particle potential obtained 
with the J04 model, cf. Fig. [9l On the other hand, they 
are largest for the most repulsive NSC89 interaction. For 
the NSC97f interaction, which is between these extremes, 
the A threshold densities are very close to those of the 
most repulsive NSC89 one, cf. Fig.[9l This stems from the 
appearance of the E~ hyperon. The effect is caused by 
the slowdown of the increase of the neutron chemical po- 
tential and is further related to the rapid increase of the 
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E~ density just after its appearance, cf. Fig.[5l Basically, 
the slowdown occurs as soon as a new hyperon appears 
because most of the energy is used for its creation. Once 
the concentration of the hyperon has reached a plateau, 
the neutron chemical potential resumes its increase until 
a further hyperon might appear. Thus, the appearance 
of the first hyperon shifts the threshold density of the 
next hyperon towards higher values. This effect explains 
why the threshold densities of the A are so similar for 
the NSC97f and NSC89 interactions. It also makes clear 
why the A threshold densities for the xEFT600 interac- 
tion are smaller than those of the NSC97a, NSC97c and 
NSC97f interactions even though their A single-particle 
potentials are almost the same, cf. Fig. [T] In the case 
of the J04 model the delay mechanism described above 
becomes very interesting. For this YN interaction the 
A and S~ hyperon appear almost at the same density. 
In this case the neutron chemical potential stagnates but 
the A and the T,~ single-particle potentials are attractive 
enough to compensate for this. 

To summarize this section strangeness appears around 
~ 2po for all YN models and parameter sets used. Note, 



that the appearance of the first hyperon, whether it is 
the I]~ or the A, cannot be altered by taking into ac- 
count YY interactions which have been neglected in this 
work. The present study in terms of the broad parameter 
ranges as well as the multitude of YN interaction models 
reveals that strangeness in the interior of neutron stars 
cannot be ignored. Similar conclusions are obtained in 
the Brueckner-Hartree-Fock theory (26| . 



V. STRUCTURE OF NEUTRON STARS 

In this section we analyze the effect of the EoS in- 
cluding hyperons on neutron stars. We focus on non- 
rotating stars, ignoring any changes, caused by the ro- 
tation. For a given EoS, the mass-radius relation of a 
NS can be determined by solving the familiar Tolman- 
Oppenheimer-Volkoff equation (TOV) To describe 

the outer crust and atmosphere of the star i.e., the re- 
gion of subnuclear matter densities for very small baryon 
densities below pB < 0.001 fm~^, we have used the 
EoS of Baym, Pethick, and Sutherland (2^, which re- 
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lies on properties of heavy nuclei. For densities between 
0.001 fm"^ < Pb < 0.08 fm"^,i.e. for the inner crust, we 
have used the EoS of Negele and Vautherin [l^ , who have 
performed Hartree-Fock calculations of the nuclear crust 
composition. Details on crust properties can be found 
e.g. in [sO, while recent state-of-the-art approaches 
are discussed in (3^ . 

In Fig. [in] the mass-radius relation of a NS for a soft 
EoS (left panel) and for a stiff EoS (right panel) is shown. 
The symmetry energy at = 32 MeV is kept fixed in both 
calculations and the resulting mass-radius relation with- 
out any strangeness is also added for comparison. 

As can be seen from the figure the appearance of hyper- 
ons reduces the NS mass drastically compared with the 
pure NN case. Even for larger values of the incompress- 
ibility, i.e. Kq = 300 MeV, the maximum mass, obtained 
with all YN interactions used, is still below the largest 
precisely known and measured NS masses 1.44 Mq of the 
Hulse- Taylor binary pulsar. This is not an unusual result 
and is also seen in other, related works such as e.g. (8l.l33l- 
[35| . In general, any inclusion of further degrees of free- 
dom to the nucleons will reduce the NS mass. 

Due to large uncertainties in the high density behavior 
of the symmetry energy we investigate its infiuence on 
the NS masses as follows: The parameter 7 in Eq. ^ 
determines the symmetry energy changes, i.e., how asy- 
stiff or asy-soft the EoS is. The value of 7 = 0.6, used so 
far, represents a asy-soft system [s^. Above saturation 
density we change the value of this parameter and use 7 = 
1 while keeping all other parameters fixed representing a 
asy-stiff system. The effect of the 7 modification is shown 
in Fig. [ni 

Only small changes in the mass-radius relation are ob- 
tained when hyperons are included. In most cases the 
mass difference is O.IMq. An increasing 7 leads to an 
increase in the proton concentration which in turn leads 
to an increase in the concentration. This largely can- 
cels any energy gain from an increased symmetry energy 
term that could possibly increase the mass of a neutron 
star. The biggest mass change is seen for Kq = 200 MeV 
with the xEFT600 interaction since this model does not 
contain any E~ hyperons. However, even in this case 
the difference is below O.2M0 which leads to a maximum 
mass below 1.2M0. A more noticeable change can be 
observed for the radius. A radius shift of 1 km towards 
larger radii is seen for all curves. Furthermore, if we in 
addition vary at at saturation density between 28 and 36 
MeV, a less pronounced effect on the NS mass is visible 
as compared to the variation of 7. 

To complete our investigation of hyperon effects on the 
neutron star EoS we also need to consider modifications 
induced by the YY and 'E.N interactions. In order to eval- 
uate the effects of the S = —2 sector we have to construct 
YY and SiV V^ow k interactions, along the lines discussed 
earlier for the S* = — 1 sector. Obviously, since there are 
more particles to consider, the situation complicates con- 
siderably from a numerical point of view. Unfortunately, 
unlike the S — —1 sector where we had several different 



"bare" interactions from which we constructed the Viow k 
potentials, in the S = —2 sector we only have one, namely 
the NSC97 interaction. We also note that the inclusion 
of the 5 = — 2 sector will not influence the appearance of 
the first hyperon. Hence we have neglected the S = —2 
up to this point. 

Fig. [T^ shows the mass-radius relation of a NS with dif- 
ferent yiV, YY and SiV interactions is presented. The 
effect of the inclusion of the S — —2 sector is rather 
marginal. As is visible in Fig. [T^ the maximum masses 
are lower than in the previous cases which is reasonable 
since a further degree of freedom, the S particle, is added. 
However we note that the NSC97 YY interaction is at- 
tractive which is the reason for the decrease of the allowed 
maximum NS masses, see e.g. Q. A repulsive YY inter- 
action would have an opposite effect as shown in js^l- 

It is also interesting to point out that for the xEFT600 
interaction, where the hyperons appears late, if at all, 
it is the S hyperons which takes its pace and influences 
the maximum mass of a NS [13, [3 • For all other models 
in which the S~ appears earlier than the S their influence 
on the NS mass is marginal. 



VI. SUMMARY AND CONCLUSIONS 

The main intention of the present work was to study 
the consequences of available hyperon-nucleon interac- 
tions on the composition of neutrons stars and the max- 
imum masses. The analysis was performed in the frame- 
work of the renormalisation group improved Viow k in- 
teraction deduced from available bare potentials. Since 
the experimental data base is very limited, these poten- 
tials are not well constrained, in contrast to the nucleon- 
nucleon case. 

We have determined the threshold densities for the 
appearence of hyperons in /3-equilibrated neutron star 
matter to lowest order in a loop expansion. To explore 
the sensitivity to available YN potentials we have con- 
structed single-particle potentials for the A and S~ hy- 
peron in lowest order and deduced the energy per parti- 
cle. We have replaced the pure nucleonic contribution to 
the energy per particle by an analytic parameterization. 
This replacement enables us to vary the incompressibility 
Kq and symmetry energy at of the purely nucleonic EoS 
and to investigate their influence on the threshold den- 
sities of hyperons. The composition of /3-stable matter 
has been determined by the requirement of charge neu- 
trality and /3-equilibrium. The corresponding threshold 
densities for various values of A'o and at were evaluated. 
The most important conclusion is that a more attrac- 
tive single-particle potential will decrease the chemical 
potential and thus decrease the threshold density of the 
corresponding hyperon. 

We have found that, irrespective of the YN interac- 
tions, incompressibility and symmetry parameter used, 
hyperons will appear in dense neutron star matter at 
densities around 2po- This inevitably leads to a signifi- 
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FIG. 10: Mass-radius relation of a neutron star for symmetry energy at = 32 MeV and different YN interactions. For 
comparison the mass-radius curve obtained for tlie pure A'^A'' interaction is also shown. Left panel: soft EoS, right panel: stiff 
EoS. 
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FIG. 11: Mass-radius relation of a neutron star for symmetry energy at — 32 MeV, 7 = 1 for p > po and different YN 
interactions. For comparison the mass-radius curve obtained for the pure A^'A'^ interaction is also shown. Left panel: soft EoS, 
right panel: stiff EoS. 



cant softening of the EoS which in turn results in smaller 
maximum masses of a neutron star compared to a purely 
nucleonic EoS. Notably, the predicted maximum masses 
are well below the observed value of 1.4 Mp) , an outcome 
also known from other works, e.g. [sl. [ssHSq . This poses 
a serious problem. 

The softening of the EoS due to hyperons cannot be 
circumvented by stiffening the nucleonic EoS, i.e., by in- 
creasing Kq, since this will cause hyperons to appear ear- 



lier. Changing the high-density behavior of the symmetry 
energy dependence or including the S = —2 sector does 
not alter this conclusion either. For more details about 
the S = —2 sector, in particular S hyperons in dense 
baryonic matter see e.g. [soMi^ . This can only mean, 
that correlations beyond the one-loop level could be im- 
portant to stiffen the hyperon contributions to the EoS. 
This, however, is not sufficient as Brueckner-Hartree- 
Fock calculations indicate H^. As has been known 
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FIG. 12: Mass-radius relation of a neutron star for symmetry energy at = 32 MeV and different YN, YY and HA'' interactions. 
For comparison the mass-radius curve obtained for the pure A''A^ interaction is also shown. Left panel: soft EoS, right panel: 
stiff EoS. 



for a long time from non-relativistic nuclear many-body 
theory, three-body interactions are crucial to yield a stiff 
nucleonic EoS. Repulsive three-body forces may also play 
a role in the hyperon sector. There is, however, little em- 
pirical information available at present. There is also the 
possibility of an early onset of quark-hadron transition 
to cold quark matter. This might also stiffen the high- 
density equation of state. 

Clearly, our results pose significant restrictions on any 
reasonable equations of state employed in the study of 
neutron star matter. With the prediction of a low on- 
set of hyperon appearance it becomes mandatory to seri- 
ously consider strangeness with respect to neutron stars. 
Eventhough our predictions for the maximum masses of 
neutron stars are too low, the treatment of hyperons in 
neutron stars is necessary and any approach to dense 



matter must address this issue. 
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